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Abstract 



The electromagnetic decay and final annihilation of magnetic monopole- 
antimonopole pairs formed in the early universe has been proposed as a pos- 
sible mechanism to produce the highest energy cosmic rays. We show that for 
a monopole abundance saturating the Parker limit, the density of magnetic 
monopolonium formed is many orders of magnitude less than that required to 
explain the observed cosmic ray flux. We then propose a different scenario in 
which the monopoles and antimonopoles are connected by strings formed at a 
low energy phase transition (~ 100 GeV). The bound states decay by gravita- 
tional radiation, with lifetimes comparable with the age of the universe. This 
mechanism avoids the problems of the standard monopolonium scenario, since 
the binding of monopoles and antimonopoles is perfectly efficient. 
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I. INTRODUCTION 



The observation of ultra-high-energy cosmic rays (UHECR) with energies above 10 11 GeV 
]l||J] poses a serious challenge to the particle acceleration mechanisms so far proposed. This 
fact has motivated the search for non-acceleration models, in which the high energy cosmic 
rays are produced by the decay of a very heavy particle. Topological defects are attractive 
candidates for this scenario. Due to their topological stability these objects can retain their 
energy for very long times and release quanta of their constituents, typically with GUT scale 
masses, which in turn decay to produce the UHECR. 

Various topological defect models and mechanisms have been studied by numerous au- 
thors ||. In this paper we investigate two different scenarios involving the annihilation of 
monopole-antimonopole pairs. We first discuss standard magnetic monopole pair annihi- 
lation [§],[5[], paying particular attention to the kinetics of monopolonium formation. We 
find that, due to the inefficiency of the pairing process, the density of monopolonium states 
formed is many orders of magnitude less than the value required to explain the UHECR 
events. 

We then present a different scenario in which very massive monopoles (m ~ 10 14 GeV) 
are bound by a light string formed at approximately 100 GeV. These monopoles do not 
have the usual magnetic charge, or in fact any unconfined flux. Gravitational radiation is 
the only significant energy-loss mechanism for the bound systems.[] Their lifetimes can then 
be comparable with the age of the universe, and their final annihilation will then contribute 
to the high energy end of the cosmic ray spectrum. 

II. REQUIRED MONOPOLONIUM ABUNDANCE. 

What density of decaying monopolonium states is required to produce the observed 
cosmic rays? The monopolonium will behave as a cold dark matter (CDM) component and 
will cluster in the galactic halo, producing a high energy spectrum of cosmic rays without the 
Greisen-Zatsepin-Kuzmin (GZK) cutoff 0,f|. Since the observational data does not seem to 
show any such cutoff, this is an advantage of such topological defect models ||[IU . 



For a given monopole mass, we can set the lifetime of the monopolonium at least equal 
to the age of the universe, and obtain the required density of monopolonium in the halo by 
normalizing the flux to the observed high energy spectrum ||. The required number density 
decreases with the monopole mass, so as a lower limit we can take the required density 
corresponding to m M = 10 17 GeV 0, 

iv^(To) >6x nr 27 cm- 3 . (1) 

Since the different components of the CDM cluster in the same way we can use this halo 
density to get the mean density in the universe, by computing, 



1 Such systems were studied in a different context by Martin and Vilenkin ||. 
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For VLcDhih 2 = 0.2, pQ DM = 0.3 GeV cm -3 , and p CI = 10 4 h 2 eV cm -3 , we get 

iWT )>10- 32 cm- 3 . (3) 

We will work with a comoving monopolonium density T = N M ^ / s where s is the entropy 
density, currently s~3x 10 3 cm -3 , so that we require 

T > 10~ 35 (4) 

to explain the observed UHECR. 

III. MAGNETIC MONOPOLE STATES 
A. Introduction 

Monopolonium states are expected to have been formed by radiative capture if there 
was a non-zero density of free monopoles in the early universe. They will typically be 
bound in an orbit with a large quantum number, so we can treat them as classical objects 
emitting electromagnetic radiation as they spiral down to deeper and deeper orbits, until 
they annihilate in a final burst of very high energy particles. 

The electromagnetic decay of monopolonium was analyzed by Hill [[| using the dipole 
radiation formula. The rate of energy loss is^ 

dt 3 g 2 M m 2 M ' 

where gu is the magnetic charge. From this expression, the lifetime of monopolonium with 
radius r and binding energy E = g 2 M /2r is [|] 

°9m 



For mjif = 10 16 GeV, qm = V(2e) ~ V%i, and an initial radius of r = 10 _9 cm, this gives 
te ~ 10 18 sec, comparable to the age of the universe. 

Bhattacharjee and Sigl used a thermodynamic equilibrium approximation to estimate 
the monopolonium density and argued that the late annihilation of very massive magnetic 
monopoles could explain the UHECR events observed. Here we recalculate the density of 
monopolonium states, taking into account the kinematics of formation and the frictional 
energy loss of monopolonium formed at early times. 



2 Here and throughout we use units where h = c = ks = 1. 
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B. Friction 



Before electron-positron annihilation, monopoles interact with a background of rela- 
tivistic charged particles. These interactions produce a force which, for a non-relativistic 
monopole is given by []TT] 

F = -s N < T v L T < 7 > 

where N c is the number of species of charged particles, v the velocity of the monopole with 
respect to the background gas of charged particles and b the impact parameter of the incident 
particles. Since we are interested in the friction that a monopole feels in a bound state orbit 
of monopolonium, we will not consider the interaction of charged particles with impact 
parameter greater than the radius of the monopolonium, so 6 max ~ g 2 M E~ x . Initially, the 
monopoles are bound with energy E ~ T, so 6 max ~ g 2 M T~ l . Equation (|7|) is derived using 
the approximation that each charged particle is only slightly deflected. This approximation 
breaks down for impact parameters that are too small, so we should cut off our integration 
at Q b min m T~ l . Using N c = 2 and g 2 M « 34, we get 

FW1.22T 2 !) (8) 

so the energy loss rate due to interactions with charged particles in the background is 

— « 1.22 TV. (9) 

dt v ; 

Taking the system to be bound in a circular orbit, we have 

m M v 2 ~ E (10) 

so we can write 

(11) 

at ttim 

The time scale for this process is 

TF dE/dt ~ 1.22 T 2 1 j 

If we compare it with the Hubble time, 



/ 90 

r " = \hr^ m pi T ~ 2 ~ 0.184m pi T" 2 , (13) 

where m p ; is the Planck mass, and g* is the number of effectively massless degrees of freedom, 
g* = 10.75, we get 

^0.15^«1. (14) 
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Thus, we see that the damping of the monopolonium energy due to friction is very effective 
in this regime, and the monopoles spiral down very quickly. 

When the distance between monopoles becomes small as compared to T _1 , the effect 
of friction is reduced and Eq. (|7|) is no longer accurate. However, even for T = 1 MeV, 
the radius has been reduced about two orders of magnitude to r ~ 2 x 1CT 11 cm, and the 
electromagnetic lifetime has been reduced by about six orders of magnitude. Thus only 
monopolonium states formed after electron-positive annihilation can live to decay in the 
present era. 

After electron-positron annihilation the number of charged particles in the thermal back- 
ground has decreased by a factor ~ 1CT 9 so t f /t h ^ 1 and the monopolonium is little 
affected by friction. 



C. Formation rate 

We can obtain an upper limit for the monopolonium density by solving the Boltzmann 
equation, 

dN M M 2 

df = {vbV)n M -3HN M m, (15) 

where tim denotes the free monopole density, N M] ^ the monopolonium density, H the Hubble 
constant, and {a b v) the average product of the binding cross section times the thermal 
velocity of the monopoles. 

With the comoving monopole density 7 = hm / s, we can rewrite the equation above as 

— = (a b v)^n M = (cr b v)rfs. (16) 

Using the approximation for the classical radiative capture cross section of monopoles with 
thermal velocities given by |12| , 



M *^T^A-T) ' (17) 



and with 

2tt 2 



9*sT\ (18) 



45 

where g*s is the number of degrees of freedom contributing to the entropy, we get 

dT vr^of^/m.A 9 / 10 



Tt ~ 45 ml VT) S * ST • (19) 



Since we are interested in the evolution of the monopolonium density after electron-positron 
annihilation, we will take a constant value g*s ~ 3.91 to get 

* K426 ^(#) 1 (20) 
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As we will see, only a tiny fraction of the monopoles will ever be bound, so we can 
consider the comoving number of monopoles 7 to be constant. To integrate Eq (|20"D , we will 
make a change of variable 



/ 90 

t = "V T ~ 2 ~ 0.1647VT- 2 , (21) 



appropriate to times after electron-positron annihilation, to get 

dT m 2 M \ T J 



dv _ 1M <!Mf(^\ (22) 



r^m^^pY'V (23) 



and thus 



\m M J \m M 

We now take T{ ~ 1 MeV and g 2 M = 34, and note that to produce the observed UHECR, we 
must have itim > 10 11 GeV, so that for a fixed monopole comoving density 7, we have the 
bound, 

T f < 4 x 10 V . (24) 



D. Monopole density bound 



The formation of magnetic monopoles via the Kibble mechanism [13] is inevitable in all 
GUT models of the early universe, and annihilation mechanisms are not efficient in a rapidly 
expanding background fl4| , [T5|| , so that the typical initial density of monopoles produced at a 
GUT phase transition will very soon dominate the energy density of the universe. The most 
attractive solution for this problem is the inflationary scenario |16[]. In standard inflation, 
the exponential expansion of the universe reduces the monopole density to a completely 
negligible value. However, it is possible for new monopoles to be formed at the end of 
inflation [jT^-[21"||. The exact relic abundance of monopoles created in this period is very 



model dependent, but its value is constrained by the Parker limit [p^] : To prevent the 
acceleration of monopoles from eliminating the galactic magnetic field, the monopole flux 
into the galaxy must be limited by 

F < 10~ 16 cm -2 s _1 sr -1 . (25) 

Assuming a monopole velocity with respect to the galaxy of ~ 10~ 3 c, we can translate this 
bound into a limit on the monopole density, 

n M < 10~ 23 cm" 3 , (26) 

and thus 7 < 10 _26 .Then, from Eq. ( p4[) we have 

T f < 10~ 45 . (27) 
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Since this conflicts with Eq. (£|) by 10 orders of magnitude, we conclude that primordial 
bound states of magnetic monopoles cannot explain the UHECR. 

We note that we have used several approximations which overstate the possible value of 
Ty: First, we have considered the total classical radiative capture cross section. This takes 
into account not only the monopolonium formed with the right energy to decay at present, 
but all the possible binding energies, clearly overestimating the value of Tf. Second, it has 
been argued that the classical cross section given in Eq (O) overestimates its real value due 



to photon discreteness effects fll2f . Finally, some of the monopolonium will have decayed 
before the present time, reducing the value of T. All of these effects make the conflict above 
more serious. 

IV. MONOPOLES CONNECTED BY STRINGS. 

We present now a different scenario for the formation and annihilation of monopole- 
antimonopole bound states. The main problem in explaining the UHECR by the conven- 
tional magnetic monopolonium system is the inefficiency of the binding mechanism. This 
can be solved if we assume that all the monopoles get connected by strings in a later phase 
transition. Since the U(l) symmetry of the monopoles would be broken by the second phase 
transition, this U(l) must be a field other than the usual electromagnetism.^ We further- 
more assume that these monopoles will not have any other unconfined charge, so that they 
will feel almost no frictional force moving in a background of particles. 

We take the co moving density of bound monopole systems T to be constant. With a 
monopole mass of 10 14 GeV the calculation of Sec. [Tl| gives T ~ 10~ 33 , and with all monopoles 
bound, 7 = 2r. The proper density at the time of string formation is then 

n M (T s ) =1 s = ^g^sTh ~ 10~ 32 T 3 . (28) 
45 

We can then compute the mean separation between monopoles at the time the string is 
formed, 

U ~ [n M (T s )]- 1/3 . (29) 
If we take T ~ 100 GeV, we obtain 

Li~10- 6 cm, (30) 

which is much smaller than the horizon distance, du ~ 3 cm at T ~ 100 GeV. We will 
assume that there are no light (m ~ T s or less) particles that are charged under the string 
flux. This means that there will be no charged particles that interact with the monopoles and 
cause the system to lose energy, so that gravitational radiation will be the only energy loss 



3 This is different from the Langacker-Pi scenario [ 23 1 , where electromagnetism is broken and then 



restored at a lower temperature, and monopoles do feel large frictional forces. 
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mechanism. When the strings are formed they may have excitations on scales smaller than 
the distance between monopoles, but these will be quickly smoothed out by gravitational 
radiation, leaving a straight string. The energy stored in the string is then /iLi, where 
[i ~ T s 2 is the energy per unit length of the string. This is smaller than the monopole mass 
by the ratio 

„ io- 2 

so the monopoles will move non-relativistically. 

In order to estimate the radiation rate we can assume that the monopoles are moving 
in straight lines. In fact, at the time of string formation the monopoles will have thermal 
velocities, so that in general the system will be formed with some non-zero angular mo- 
mentum. However, in general this will be small compared to the linear motion due to the 
string tension, so we will ignore it, except to note that the monopoles will pass by without 
collision. The half oscillation of one monopole is parameterized by 

x(t) = (2aLf/H - \at 2 (32) 

with a = n/wriM and < t < (8-L/a) 1 / 2 . Using the quadrupole approximation,^ the rate of 
energy loss of the system is 

dt 45 \m M J 
Since \iL is the energy in the string, we can integrate this equation to obtain 

L = Ue- 11 ^ 

with 

45 m M 45 m^rriM 
Tg ~ 288CV ~ 288 T s 4 

The lifetime of the state will thus be r g ln(Lj/r M ), where r M ~ mj} is the radius of the 
monopole core. For T ~ 100 GeV and tum ~ 10 14 GeV, Eq. ([35]) gives r g ~ 10 17 sec, 
comparable with the age of the universe. 

This suggests that the bound system formed by a monopole-antimonopole pair connected 
by a string can slowly decay gravitationally, and release the energy stored in the monopole 
in a final annihilation when the two monopole cores become close enough. 



4 The fully relativistic situation was considered in Q. 
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(31) 



(33) 



(34) 



(35) 



V. CONCLUSIONS 



We have shown that is not possible to construct a consistent model for the origin of the 
UHECR based on the electromagnetic decay and final annihilation of magnetic monopole- 
antimonopole bound states formed in the early universe. We have obtained an upper limit 
for the monopolonium density today, taking into account its enhancement in the galactic 
halo and the maximum average free monopole density consistent with the Parker limit. Due 
to the small radiative capture cross section for the monopoles and the rapid expansion of 
the universe, the maximum density of monopolonium is many orders of magnitude below 
the concentration required to explain the highest energy cosmic ray events. 

We then proposed a different scenario in which the monopoles are connected by strings 
that form at a relatively low energy. This mechanism solves the problem of the inefficiency 
of the binding process, since every monopole will be attached to an antimonopole at the 
other end of the string. Due to the confinement of the monopole flux inside of the string , 
the main source of energy lose for these bound systems will be gravitational radiation. If 
we assume a monopole mass of 10 14 GeV and a string energy scale of the order of 100 GeV, 
the lifetime of the bound states would be comparable with the age of the universe, making 
them a possible candidate for the origin of the ultra-high-energy cosmic rays. 
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